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ABSTRACT 

To detect Earth-like planets in the visible with a coronagraphic telescope, two major noise sources have to be 
overcome: the photon noise of the diffracted star light, and the speckle noise due to the star light scattered by 
instrumental defects. Coronagraphs tackle only the photon noise contribution. In order to decrease the speckle 
noise below the planet level, an active control of the wave front is required. We have developed analytical methods 
to measure and correct the speckle noise behind a coronagraph with a deformable mirror. In this paper, we 
summarize these methods, present numerical simulations, and discuss preliminary experimental results obtained 
with the High-Contrast Imaging Testbed at NASA's Jet Propulsion Laboratory. 
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1. INTRODUCTION 

In coronagraphs like any other optical system, aberrations give rise to scattered light in the form of speckles 
in the focal plane, thus preventing very high-contrast imaging. A wave front sensing and control system can 
remedy this, and make the detection of faint objects possible.^ Such a system is absolutely critical to the success 
of NASA's Terrestrial Planet Finder Coronagraph (TPF-C)^ that aims not only at the detection, but also at the 
spectroscopy of terrestrial planets in the visible range. This is a very challenging task as the contrast* between a 
terrestrial planet and its Sun-likc star amounts to ~ 10^" in the visible. With that goal in mind, we have devised 
a theory'^ to measure and correct wave front aberrations, so as to clear out of speckles an area of the focal plane, 
referred to as dark hole (DH). This area is then suitable for very faint object detection. 

In the following, we first summarized the theory {^^, then we present simulations {^^, and finally we discuss 
preliminary experimental results (2)1. 

2. THEORY SUMMARY 

The wave front sensing and control system we advocate consists of 

1. The sensor: the science camera itself since a separate sensing channel would have its own aberrations, 
impossible to disentangle from the main optical train's aberrations f ii2.1|) : 

2. The actuator: a deformable mirror (DM) with NxN individual actuators f i)2.2|l : 
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3. The command control system: a computer converting speckle field measurements into DM commands 




Let us now examine briefly these three elements. 

2.1. Speckle field measurement 

In a regime of small aberrations, the electric field in the image plane after the coronagraph is the sum of two 
terms: the original speckle field, $, and the DM-controUed electric field, Xf. In principle, three images with 
different DM settings bring enough information to measure the complex amplitude of the speckle field. Indeed, 
the intensity /„ (n = 0, 1, 2) in any given pixel of image n reads 

(1) 

[ h = \^ + ^0 + 5^2? ■ 

System ^ has for solution 

with the caveat that the denominator should not be zero. This tells us that the successive modifications in the 
DM-controUed electric field, S'i/i and ^^'2, must lead to an intensity change in any given pixel big enough so 
that information could be acquired for that pixel during the measurement process. 

2.2. Deformable mirror model 

Assuming linear superposition and introducing actuator influence functions^, the DM-controUed electric field can 
be expanded as 

fc=0 1=0 

where aki are the actuator strokes and F^i the complex images on the detector of the influence functions. For 
the simulations and the experiments, we used a 32 x 32 DM with a 1-mm actuator spacing. More information 
about the DM can be found in Ref. 4. 

2.3. Speckle nulling 

Field nulling In this approach, we try to cancel out the speckle fleld in every detector pixel. The basic 
equation is simply $ + = 0, i.e. by making use of Eq. Q, 

N-lN-l 



auiFu = (4) 



fc=0 (=0 



Equation Q is a linear system in the actuator strokes that can be solved either by singular value decomposition,^ 
or very efficiently with an inverse Fast Fourier Transform (FFT), provided Fki are phase-shifted copies of one 
another. Although this last hypothesis does not rigorously hold behind a coronagraph, simulations show that 
the solution by inverse FFT could yield satisfactory results for a first round of tests. 

^An influence function describes the DM surface deformation in response to the actuation of a given actuator. 



Energy minimization In this approach, we seek to minimize tlie total energy of the speckles in the DH. This 
energy is defined as 

S= [f + (5) 



Equation @ is minimized when the energy gradient with respect to the actuator strokes is zero, i.e. 



n— m— 



Equation (jSJ is again a linear system in the actuator strokes, but with smaller matrices than in equation |0J 
because matrix sizes are governed here solely by the number of actuators, and no longer by the number of 
detector pixels. Simulations show that speckle energy minimization is more powerful and flexible than speckle 
field nulling. We intend to demonstrate experimentally this second approach in the near future. 



3. SIMULATIONS 

In this section, we present monochromatic simulations using the field nulling equation (Eq. ^ solved with an 
inverse EFT (Eigs.^121). We consider a Lyot-type coronagraph with a band- limited image-plane mask.^ Because 
the DM features 32 x 32 actuators on a square grid, the sampling theorem imposes that the maximum size for 
the DH is a 32-^ X 32-^ square, where A is the wavelength and D the entrance pupil's diameter. When both 
phase and amplitude defects are present, a single DM creates a DH limited to half the maximum size. The 
simulations in Figs. are meant to illustrate the theory, not to be realistic. They assume a Gaussian white 
noise for phase defects and amplitude defects, but no photon nor detector noise. The contrast in any given pixel 
is defined as the ratio of the intensity in that particular pixel to the peak intensity measured in the image when 
the coronagraph image-plane mask is removed. In the next section, we discuss preliminary experimental results 
and compare them with more realistic simulations. 



4. EXPERIMENTS 

In this section, we show monochromatic experiments at A = 785 nni with the High Contrast Imaging Testbed^ 
at JPL (Eigs.|31El). For these, we have used the field nulling equation (Eq. Q solved with an inverse EFT. In the 
initial image, the median and the mean in the DH are 73 and 146 ADU, respectively, corresponding to contrasts 
of 1.4 X 10~^ and 6.8 x 10"^, respectively. In the last image, the median and the mean in the DH are 72 and 
139 ADU, respectively, corresponding to contrasts of 1.3 x 10"^ and 6.7 x 10^^, respectively. Unfortunately, 
the contrast remains almost unchanged in this first attempt. Nevertheless, it makes sense at this stage to try to 
iterate the measurement and correction process because some modeling approximations can be regarded as phase 
and amplitude errors. Indeed, after five iterations, the median and the mean in the DH are 51 and 112 ADU, 
respectively, corresponding to contrasts of 1.0 x 10~^ and 4.8 x 10^^, respectively. Simulations with a realistic 
noise level show that more iterations could have led to a median DH of ~ 5 ADU, about ten times better than 
what was obtained. 



5. CONCLUSION 

Although our first experimental results are modest, they prove that our approach is sound and should be pursued. 
Future work will include an experimental demonstration of the energy minimization approach in monochromatic 
light, as well as a generalization of the theory to polychromatic light with corresponding experiments. Our wave 
front sensing and control theory is general and can be applied to other types of coronagraph, as demonstrated 
in Ref. 8 for shaped-pupil coronagraphs. 

Note that these experiments are by no means representative of the best result achieved to date with the 
HCIT. For these, we refer the reader to Ref. 7 and to the paper by Trauger et al. in these proceedings. 




Figure 1. Simulated series of deformable mirror (DM) settings during the speckle nulling process. The first image on the 
left displays the initial DM setting (dmO), whereas other images (dml, dm2, & dm3) show changes with respect to that 
initial setting. Note that the DM is not initially flat because the bulk of low-order aberrations has already been taken 
out. In the second and third steps (dmf & dm2), the speckle field is modified so as to measure its complex amplitude. 
In the fourth step (dm3), the shape imposed on the DM creates the dark hole. The DM features 32x32 actuators, but 
actuators outside the entrance pupil (a 30-mm diaphragm on top of the DM) are not commanded. Actuator strokes are 
color-coded in nm. 




Figure 2. Simulated series of detector images during the speckle nulling process. Every image displays the log of the 
contrast as defined in text. Superimposed is a dotted grid scaled in elements of resolution from ~32-^ to -1-32-^. The top- 
left image shows the initial state: the speckle field covers the whole image except a vertical streak where the coronagraphic 
mask is opaque. In the second (top-right) and third (bottom-left) images, the speckle field is modified so as to measure 
its complex amplitude. In the fourth (bottom-right) image, the dark hole is apparent. Only the right side is corrected for 
both amplitude and phase aberrations. 




Figure 3. Same as Fig.Qbut for a real experiment. 




Figure 4. Same as Fig.|21but for a real experiment. The white rectangle ([4-^, 12-^] x [—12-^, 12-^]) delimits the region 
where the algorithm is trying to create a dark hole. The initial image already contains an annular DH that was created 
by a different method. 
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